Search for Multiple Step Integer Quantum Hall Transitions 
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Recent experiments in the integer quantum Hall regime seem to find direct transitions from a 
quantum Hall state with Hall conductance a xy — ne 2 /h with integer n > 1, to an insulating state 
in weak magnetic fields. We study this issue using a variation of the tight-binding lattice model 
for non-interacting electrons. Although quantum Hall transitions with change in Hall conductance 
ne' 2 jh with n > 1 do exist in our model for special tuning of parameters, they generically split 
into quantum Hall transitions with the Hall conductance changing by e 2 jh at each transition. This 
suggests that a generic multiple step quantum Hall transition is incompatible with a non-interacting 
electron picture. 



PACS numbers: 73.43. Cd, 73.43. Nq, 71.30.+h 

For non-interacting two-dimensional electron systems 
in the presence of a perpendicular applied magnetic field 
(B), Kivelson, Lee and Zhang [Q proposed a global phase 
diagram, which predicted that in a spin-polarized system, 
at a quantum Hall transition, the Hall conductance (o~ xy ) 
can only change by e 2 jh as filling factor (n) changes in 
the integer quantum Hall effect (IQHE) regime. Stimu- 
lated by Kivelson et a/.'s proposal, a number of experi- 
ments reported the observations of a xy = — » 2e 2 /h — > 
transition [|] ||| i which is consistent with the global phase 
diagram, assuming the Zeeman splitting in the lowest 
Landau level is unresolved. However, more recently, 
many experiments |^-|^] , some in systems with spin-split 
Landau bands, observed direct transitions from an insu- 
lating state at low fields to a xy — ne 2 /h IQHE states with 
n upto 6. Henceforth in this letter, we dub these tran- 
sitions from the insulating phase to n > 1 IQHE states 
multiple step quantum Hall transitions (MSQHTs). 

In a series of papers |9j-fh2]], Sheng et al. claimed that 
the experimental studies of the quantum Hall transitions 
in weak magnetic fields are consistent with their numeri- 
cal calculations on a tight-binding lattice model. In par- 
ticular, the MSQHTs are the consequence of the critical 
energy in the lowest Landau level floating and merging 
with critical energies in higher Landau levels How- 
ever, it remains unclear how these critical energies can 
merge, and thus give rise to an MSQHT. An alternative 
interpretation of the apparent MSQHT is a series of sin- 
gle step (Aa xy = e 2 /h) plateau transitions that are too 
close to distinguish numerically. 

The two scenarios may be resolved by studing the crit- 
ical behavior in the vicinity of the transition(s) under 
investigation. For a single, spin-split Landau level, nu- 
merical studies by various groups fll3|-[l6[ using different 
microscopic Hamiltonians are consistent with the same 
critical behavior at the quantum Hall transition (QHT). 
In all cases, the localization length diverges with an ex- 
ponent v ~ 2.3 ± 0.1 at the critical energy, E c , which 
is at the band center in the case of electron-hole sym- 
metry. By studying the localization length exponent, for 



instance, Lee and Chalker JT7| demonstrated, in a spin- 
degenerate Landau level, the existence of a pair of de- 
localization transitions in the same universality class as 
the spin-split system. 

In this letter, we report our search for MSQHTs by 
truncating the Hilbert space of the tight-binding lat- 
tice model to a single (central) magnetic subband, which 
contains a multiple of 2e 2 /h Hall conductance, and can 
therefore support a MSQHT in principle. While this 
model does not correspond directly to the experimen- 
tal situation, it allows for a much clearer interpretation 
regarding the possible existence of a multiple step tran- 
sition. We also expect the study to lead to a clearer 
understanding of the tight-binding lattice in the context 
of weak magnetic fields. 

Our results can be summarized as follows: we find that 
the system can support a MSQHT in the limit of infinite 
disorder, but in this limit the multiple step comes from 
a set of noninteracting, single step transitions with acci- 
dental degeneracy of the step position. As soon as the 
disorder is made finite, the mutiple step breaks into sin- 
gle steps, whose critical behavior is the same as that for 
isolated Landau levels. 

We consider a tight-binding model of non-interacting 
spinlcss electrons on two-dimensional square lattice, with 
nearest neighbor hopping (t), a uniform perpendicular 
magnetic field B, and on-site disorder. The Hamiltonian 
can be written as: 



' c\cj + h.c.) + ^Cicjci, 



(1) 



where cj (q) is fermionic creation (annihilation) opera- 
tors on site i, (ij) are nearest neighbor sites, t is the hop- 
ping strength, and are phase factors due to the mag- 
netic field B. Ci is a random on-site potential, with a dis- 
tribution P(e) = (27rcr 2 )" 1 / 2 exp (-e 2 /2cr 2 ). The mag- 
netic flux per unit cell is chosen to be <j)/(f>o = Y] n = 
\/{2N + 1) with integer N, where 0o = hc/e is the flux 
quantum. 

In the presence of B, the tight-binding energy spec- 
trum splits into (2N+1) magnetic subbands. Each of the 
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27V side subbands carries a Hall conductance a xy — e 2 /h, 
while the center subband carries a xy — —2Ne 2 /h. In 
contrast to other approaches |p|- |T2| , p^[ , we first truncate 
the Hamiltonian to the Hilbert space spanned by the 
eigenstates in the center subband in the absence of disor- 
der, before adding the disorder. This results in a single 
subband with o~ xy = —2Ne 2 /h for any disorder strength. 
The negative sign of a xy is not essential to the discus- 
sion here, since we are primarily interested in the scaling 
behavior near the critical point (s) of the center subband. 

We calculate, within the center subband, the Thouless 
number Q , as a measure of the diagonal conductance 
(7 xx in units of e 2 /h, for different lattice sizes L x L (L 
is the number of sites in each dimension): 



9l{E) = 



AE 



(2) 



where (|<5-E|) is the average shift in the energy level due 
to a change of boundary condition from periodic to anti- 
periodic, and AE = l/L 2 p(E) the mean energy level 
separation, p(E) being the density of states. gh{E) is 
dimensionless, and by finite-size scaling, is expected to 
be a function of L only through the ratio L/£, where £ is 
the localization length characterizing the wavefunctions 
at the energy E, which diverges at critical energies E c as 
£ ~ \E — E c \~ v . As a consequence, the square root of the 
second moment of <?l(-E), 



W{L) 



CjE-E c ) 2 g L (E)dE 



1 1/2 



(3) 



is expected to scale as W(L) ~ L~ x l v in the thermody- 
namic limit for a single critical energy. 

We first consider the case 4>/4>q — 1/3, where the cen- 



tral subband carries a xy = —2e 2 /h, and to further sim- 
plify the situation, we begin with the infinite disorder 
limit (o~/t — > oo). Equivalently, we set t = after trun- 
cating the Hilbert space to the center subband only, so 
that the energy scale is finite. Figure [l] shows Ql(E) for 
different (even) L, averaged over 250 to 2500 samples, 
depending on L. The inset shows the width W(L) of the 
curves vs. L (solid squares for even L) on a double loga- 
rithmic plot. The straight line fit implies a single critical 
energy E c — 0, and a critical exponent v = 2.2 ± 0.1, 
consistent with the value obtained for a single step quan- 
tum Hall transition in a sideband |po| or an isolated Lan- 
dau level fTg(| . The data are well fit by a gaussian form, 
9L (E) = g exp{-E 2 /2[W(L)} 2 ), with g Q ~ 0.38, approx- 
imately twice the value go ~ 0.2 in a side subband |po| , pT| . 

To understand this result, we examined the structure 
of the projected Hamiltonian. Upon application of the 
projection operator, the on-site random potential be- 
comes non-local, giving rise to (short-range) hopping. 
In the t = limit, there is no direct hopping, and the 
hopping generated by the random on-site potential due 



to projection, leads to a bipartite Hamiltonian. Thus, 
for even-sized lattices, the system splits into two dis- 
connected sublattices interpenetrating each other, with 
identical statistical properties, since hopping from one 
sublattice to the other is not generated. This explaines 
why we see a single QHT at E = 0, with the same value 
of v as a single-step QHT, but twice the Thouless con- 
ductance go - we simply have two independent transitions 
at the same energy. 
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FIG. 1. Thouless number gh(E) of the center subband 
(N = 1) with infinite disorder (a = 1 and t = 0) for lattices 
with size L = 18, 30 and 42. The inset shows the width W(L), 
defined by Eq. (g), of corresponding gh{E), on log- log scale. 
Note the difference between even L (solid squares) and odd 
L (open squares); as explained in the text, when shifted by 
multiplying L by a factor of v2, results for odd-sized lattices 
(open circles) fall onto the dashed line fitted to even-sized 
lattices (solid squares). 

The above result is confirmed by the observation that 
for odd L, where hopping from one sublattice to the other 
is allowed by the boundary conditions, the width W(L) 
are quite different from those for even L (see solid squares 
in the inset of Fig. [l]. This seemingly strange result can be 
understood by rearranging one sublattice by translating 
its members appropriately by L lattice constants along 
one of the sides of the original lattice. One thus "unfolds" 
an odd-sized lattice to a new square lattice rotated by n/4 
and enlarged by y/2 in linear size. By replotting W(L) 
of odd L as a function of y2L, which is equivalent to 
translating the original points to the right on the log- 
log plot (Fig. H inset), we found that these points fall 
onto the same straight line obtained from those of even- 
sized lattices. Therefore, the Hamiltonian in an odd- 
sized lattice describes a single quantum Hall transition 
on a "folded" lattice. Study of statistical properties of 
energy spectrum have also confirmed that the odd-even 
anomaly is the consequence of whether the Hamiltonian 
is bipartite |p0[ . 

We next consider the more generic case, with finite dis- 
order. If, in the thermodynamic limit, there remains a 
single critical energy at the band center, one would ex- 
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pect curves for W(L) versus L to be straight lines on a 
log-log plot for large sizes, with a universal slope, related 

to the localization length critical exponent for a double spin-orbit coupling fl 
step QHT. Instead, we find (see Fig. [2]) a series of straight 
lines with a slope that depends on the disorder (a/t); fur- 
thermore, the peak value, #l(0) is dependent on a/t as 
well. (Note that the odd-even effect gets smaller as a/t 
decreases). Non- universal v for other disorder distribu- 
tions P(e) have been observed as well 



disorder case, leading to a behavior similar to that ob- 
tained in double-layer systems (2^] or for electrons with 
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FIG. 2. Width of Thouless number curves, W(L), as a 
function of lattice size L for a/t = 1, 2, 4, and 8. 

Figure |^(a) shows the full Thouless number curves for 
a/t = 5 for L = 15,21,30, averaged over 120-600 sam- 
ples; the curves show a distinct departure from the gaus- 
sian form of the scaling curve for the single-step QHT, 
with a pronounced flat-top shape. 

The apparent non-universal exponent is based on the 
assumption that the critical energy remains at E = 0; if 
instead, there were two single step QHT at ±E C in the 
thermodynamic limit, one would expect the gh{E) curves 
to be made up of two contributions. Indeed, a fit of the 
form 

g L {E) = g{L x l v \E - E c \) + g(L l ^\E + E c \), (4) 

is found to work extremely well for all \E\ > E c for all 
values of the disorder studied with a disorder dependent 
critical energy E c with a single value of v. The universal 
exponent obtained from such a fit is v = 2.3 ± 0.2, which 
is just the result for a single step QHT. Furthermore, the 
shape of the scaling curve g(L x l v \E — E c \) (see Fig. ||(b)), 
is very well fit the gaussian form obtained for a single step 
QHT. Finally, the dependence of E c on a, shown in the 
inset of Fig. |(b), is found to be E c /a = 0.78 (t/a) - 67 
for 0.03 < a/t < 1.0, consistent with E c /a — ► in the 
limit of t/a — » 0. The fit of Eq. (^) does not work as well 
between the two critical energies presumably because of 
the extremely large localization lengths, much beyond 
the sizes of systems studied, for this region. Thus the 
data implies that the accidental double step for infinite 
disorder splits into two single steps for the generic, finite 
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FIG. 3. (a) Thouless number Ql{E) of the center sub- 
band (N = 1) of lattices of linear dimensions L — 15, 21 and 
30 for a/t = 5. (b) Scaling curve of the Thouless number, 
g(L 1/u \E-E c \), with E > E c = 1.3 and v = 2.5. Dashed line 
shows the gaussian curve found for single step quantum Hall 
transitions. Inset shows that the E c , in units of a, can be fit 
to a power law (dashed line) as a function of t/a. 



We next consider the case with </>/0o = 1/5 {N = 2), 
in which the center subband has total Hall conductance 
—4e 2 /h. Since this case allows for many more possibil- 
ities, we restrict our attention to the infinite disorder 
(t = 0) case. In this limit we again obtain a Hamiltonian 
that is bipartite for even-sized lattices. This implies that 
the transition in the center subband can be either (a) 
two double step quantum Hall transitions at E = 0, or 
(b) a pair of split, single step quantum Hall transitions, 
each pair occurring at E = ±E C . Figure ^ shows gh{E) 
of even-size lattices for infinite disorder, averaged over 80 
to 1800 samples, depending on system size. 

Forcing the width W(L) to scale as a power law in L 
over the sizes studied yields an exponent v — A.2. How- 
ever, the shape of the gL(E) curves is rather flat-topped 
compared to a gaussian, and the almost doubled value 
of v also suggests that a more natural scenario is (b) 
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above. Indeed, fitting to the form given by Eq. (0) for 
\E\ > E c yields E c = 0.3, 50 = 0.4, v = 2.4, and further- 
more a gaussian shape for the scaling curve (see inset of 
Fig 4). Thus, with a single adjustable parameter, E c , 
the data are consistent with a pair of degenerate single 
step transitions, with exactly the same Thouless number 
peak, critical exponent, and scaling curve. This splitting 
of the critical energy in each sublattice argues strongly 
for single step quantum Hall transitions being the generic 
situation for non-interacting electrons. 
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FIG. 4. Thouless number ql{E) of the center subband 
(N = 2) with infinite disorder (cr = 1, t = 0) for lattices with 
i = 20, 30, 40, and 50. The inset shows the scaling collapse 
of the Thouless number g{L L ^ v \E — E c \), using the four g(E) 
curves for E > E c = 0.3 and v = 2.4. 

Our study has raised questions about the interpreta- 
tion of the results found in the full tight-binding lattice 
model |)-12|. In particular, a genuine MSQHT must be 
associated with degenerate critical energies, which is less 
likely to appear in the full model with finite amount of 
disorder and critical energies in all side subbands that 
can interact with those in the center subband. Very re- 
cently, Koschny et al. |23j have claimed to resolve two 
separate transitions with filling factor n — 2 in the full 
lattice model, and thus contradicted the direct transi- 
tion from insulator to high Hall plateaus suggested by 
Sheng et al. Nevertheless, the difficulty of unambigu- 
ously resolving this issue, as well as the obstacle that 
forbids us to obtain the scaling behavior between adja- 
cent critical energies, is consistent with the existence of 
the large localization length well beyond the sizes acces- 
sible to present numerical studies. The length scale may, 
indeed, agree with the large crossover length argued by 
Huckestein [M in the framework of the standard scaling 
theory of the IQHE. 

To summarize, we have found that for non-interacting 
electrons, even in a single band with total Hall conduc- 
tance a multiple of 2e /h, in the presence of disorder, 
a sequence of single step quantum Hall transitions ap- 
pears to be the generic situation. The only multiple step 



transition we have been able to obtain is by special tun- 
ing of parameters, for which there exist accidental de- 
generacies of critical energies along with zero coupling 
of these degenerate critical energies. In such cases, the 
multiple step transition is a trivial superposition of single 
step transitions, with the same exponents for the localiza- 
tion length (y = 2.3) and scaling curves for the Thouless 
conductance. This calls into question the interpretation 
of the apparent multiple step transitions seen in experi- 
ment j5|-|8j using a non-interacting model p|-[T2|. 

This work was supported by NSF DMR-9809483. We 
thank Kun Yang for discussions. 
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